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In hydrodynamical modeling of heavy-ion collisions the initial state spatial anisotropies translate 
into momentum anisotropies of the final state particle distributions. Thus, understanding the 
origin of the initial anisotropies and quantifying their uncertainties is important for the extraction 
of specific QCD matter properties, such as viscosity, from the experimental data. In this work 
we study the wounded nucleon approach in the Monte Carlo Glauber model framework, focusing 
especially on the uncertainties which arise from the modeling of the nucleon-nucleon interactions 
between the colliding nucleon pairs and nucleon-nucleon correlations inside the colliding nuclei. 
We compare the black disk model and a probabilistic profile function approach for the inelastic 
nucleon-nucleon interactions, and study the effects of initial state correlations using state-of-the- 
art modeling of these. 
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1. Introduction 



In ultrarelativistic heavy-ion collisions performed at the Relativistic Heavy Ion Collider (RHIC) 
and Large Hadron Collider (LHC) significant azimuthal momentum distribution anisotropies have 
been measured [jl|, |J|. These anisotropies can be explained with relativistic hydrodynamics: the ini- 
tially produced QCD-matter contains spatial anisotropies and during the hydrodynamical evolution 
these anisotropies are transferred to the momentum distributions of final state particles. 

Simulations with viscous hydrodynamics have shown that the shear viscosity of the QCD- 
matter produced in ultrarelativistic heavy-ion collisions can be estimated from the final state mo- 
mentum anisotropies. Since the origin of these anisotropies is in the initial state, uncertainties re- 
lated to the initial state must be charted before reliable estimates for the viscosity can be made. In 
this work [||] we consider two sources of uncertainties related to the Monte Carlo Glauber (MCG) 
model which is often used to initialise the hydrodynamical simulations (see e.g. [ffl], [Sj, [^]). 

The Glauber model [f7|] is usually a key element in computing the initial states for hydrody- 
namical modeling of ultrarelativistic heavy-ion collisions. Some years back, most hydrodynamical 
calculations assumed smooth initial states where the (energy or entropy) densities were assumed to 
scale with the density of binary collisions or wounded nucleons computed from the optical Glauber 
model; see, e.g. [|8j|. Now that the importance of the initial density fluctuations has been realized, 
Monte Carlo Glauber (MCG) modeling has become more frequently used. So far the black disk 
(hard-sphere) modeling of the nucleon-nucleon (AW) interactions has been the standard choice 



10), [11J] in these studies, although also more involved probabilistic ways to model the AA 



interactions have been known for a long time |2|, 13, 14]. In Ref. [||] we studied black disk and 
profile function models. 

In the MCG modeling one needs to know the positions of the initial state nucleon configura- 
tions. In most cases the nucleon positions inside the colliding nuclei are just sampled using the 
Woods-Saxon potential, neglecting nucleon-nucleon correlations [15]. However, there exist calcu- 
lations which show that high-momentum components of the nuclear wave function are generated 
by the two-body AA correlations [16, 17]. 

In Ref. H| we studied for Au-Au collision at RHIC center of mass energy y / s wv =200 GeV, 
two different uncertainties in computing the initial state asymmetries from the MCG model: one 
related to the modeling of the inelastic AA collisions between nucleons from different nuclei, 
and one related to the AA correlations in the nucleon configurations in each of the colliding nu- 
clei. In this contribution we extend the results from Ref. [||] about the anisotropy moments £„ 
with n = 1,2,3,4,5: dipole asymmetry, eccentricity, triangularity |Q, quadangular and pentagonal 
asymmetries. Results will be shown for Au+Au collisions at RHIC with an inelastic nucleon- 
nucleon cross section Onn = 42 mb. 



2. MCG framework: nucleon configurations 

The initial state of a nucleus in the MCG calculations is usually taken as a collection of parti- 
cles distributed according to a probability distribution given by the corresponding (Woods-Saxon) 
number density distribution measured in electron scattering experiments. Given the complexity 
of the nuclear many-body problem, the effects of spatial, spin and isospin dependent correlations 



2 



Initial state anisotropies uncertainties in MCG 



M. Alvioli 



among the nucleons are usually overlooked and the nucleons are positioned randomly for each of 
the simulated events. Recently, in Ref. [ |l8[ ] it was shown how such an approach can be modified 
by including initial states, which are prepared in advance, in the commonly used computer codes. 
Also, a method to produce such configurations was introduced. The method is based on the no- 
tion of a nuclear wave function \jf, which contains the nucleonic degrees of freedom and which 
is used to iteratively modify the positions of randomly distributed nucleons using the Metropolis 
method so that the final positions correspond to the probability density given by \y\ 2 . The method 
is constructed to reproduce the same nucleon number-density distribution as the usual one and, 
in addition, to reproduce the basic features of the two-nucleon density in the presence of the AW 
correlations. The model wave function is taken in the form 

A 

Y(xi,...,x A ) = Y\f i j<j>(x 1 ,...,x A ), (2.1) 

where is the uncorrected wave function and fij are correlation operators; here, xt denotes the 
position, spin and isospin projection of the i-th nucleon. The correlation operator contains a de- 
tailed spin-isospin dependence. In the most general case, this dependence includes a number of 
channels that are the same as the one appearing in modern nucleon-nucleon potentials used to 
successfully describe a variety of properties of light and medium-heavy nuclei within different ab 
initio approaches. The major merit of this approach is that the two-body densities resulting from 
configurations obtained using these correlations are clearly more realistic than the completely un- 



correlated ones Q19[]. In this paper we have used configurations generated by two-body correlations, 
including the tensor operator, and three-body clusters surrounding each of the nucleons, induced 
by full correlations. The effects of genuine three-body correlations have been discussed in Ref. [Q]. 

2.1 MCG framework: modeling the inelastic interactions 

In each simulated event, given the impact parameter of the A +A collision, the nucleon-nucleon 
interactions must be modeled. We work in the Glauber model framework [0], neglecting the effects 



of inelastic diffraction that lead to fluctuations of the strength of the AW interactions Q20| , |21 |. To 
generate the inelastic NN collisions of interest here, we use the following two different approxima- 
tions for establishing whether a collision between the nucleons i and j from different nuclei takes 
place: 

• Black disk approximation, used recently, e.g., in Ref. [Q], where one assumes the two nucleons 
to interact inelastically with a probability one if their transverse separation bij is within a radius 
defined by the inelastic NN cross section aff N , 

bjj<^; (2.2) 

• Profile function approach, where the probability of an inelastic interaction between the nucleons 
i and j is given by 

p(b iJ ) = i-\i-r(b iJ )\ 2 , (2.3) 

and where the profile function T is expressed in terms of the total and elastic NN cross sections as 
follows: 

mj) = ^e-^ 2B \ (2.4) 
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withfi = (a^) 2 /(167T<). 

The probability distribution P(bjj) can be derived in the Born approximation of the potential- 
scattering formalism [[7], 22]; details are given in Ref. [J3[. The nucleon-nucleon elastic, total 



and inelastic cross sections are given in this formalism by: 

i a tot ^2 



Onn - mnB , (2-5) 



°N*N — °NN + °NN i ( 2 -6) 

= J d 2 b u (l - |1 - r(ft ;7 )| 2 ) . (2.7) 

After the nucleon-nucleon interactions have been determined, we can calculate the initial 
asymmetries from the positions of the nucleons which had experienced at least one collision. These 
are called wounded (or participant) nucleons. 

2.2 Spatial asymmetries and their fluctuations 

We compute the spatial asymmetries from the wounded nucleon positions which are obtained 
from the MCG model as follows: 

_ (w(r)cos(w(0-y„)) 

£n — / — TTT > U-°) 

(w(r)) 



where w(r) is a weight and we choose w(r) = r 2 for n = 1,2,3,4,5, following Ref. [23] (at variance 
with Ref. [^] where we used w(r) = r 3 for n = 1 [ |26| ] and w(r) = r" for n = 2, 3). 
The orientation angle Xjfn is determined as 

1 (w(r)sm(n<b)) n 

Y„ = - arctan ) ) ( ) v » + - , (2.9) 

where the arctan must be placed in the correct quadrant. Since the experimental methods typically 
measure the root mean square (rms) of the flow coefficients, we also present the rms values of 
initial state anisotropies y(l£) 

3. Results and discussion 

In Ref. [||] we charted some of the uncertainties in the computation of the initial state ani- 
sotropies from the Monte Carlo Glauber model. A summary of the results is plotted in Fig. [j], 



now using rms values. Also harmonics n = 4, 5 [|24j, |25j], which were not shown in the original 
article, are shown here. We used two different ways of modeling the inelastic interactions between 
the colliding nucleons. The difference between these two cases gives us an estimate about the 
uncertainties related to this part of the model: in central collisions the details of the interaction 
model play a minor role, but in the peripheral collisions such details can cause uncertainties up to 
10% in the first three harmonics £i, £2 and £3. The situation is similar for the rms of £„, as shown in 
Fig. [ll where n = 1,2,3,4,5 are considered with the choice of the weight function, w(r) = r 2 ; the 
effect is at about 10% over the N part range and it is mostly due to the probabilistic profile function 
approach. We also checked in Ref. [Q] that with these two interaction models the difference in 
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Figure 1: Left panel: the root mean square of the initial anisotropies, 



defined in Eq. (2.8), as a 



function of the number of participants. Results are relative to the black disk, uncorrelated approximation 
(dotted lines) and probabilistic interaction plus correlated configurations (solid lines). Right panel: the ratio 
of the r+correlations over black disk, uncorrelated results shown in the left panel. 



the number of wounded nucleons and binary collisions remains small in central collisions, but at 
impact parameters 10-15 fm the difference can be around 10%. 

We also presented a study of the effects of AW correlations with an update of correlated config- 
urations and extended discussion as compared with the previous published papers on this subject. 
We confirmed that the inclusion of centrally correlated nucleon configurations produce the effects 
to eccentricity and its relative variance as was claimed by Ref. [^8[j. As a new result, we observed 
that the inclusion of realistically correlated configurations (two-body full correlations, three-body 
chains) seems to essentially cancel this effect and bring the results back close to the no correlations 
case. The effect is similar for dipole asymmetry and triangularity as for eccentricity. However, we 
also showed that there are still uncertainties caused by the truncation done in the nucleon configu- 
ration calculation with full correlations and we expect three-body correlations to play a role. 

In this Proceedings contributions, we added the study of the root mean square of initial state 
anisotropies for n = 1,2, 3,4,5 and presented results in Fig. [TJ with similar sensitivity to the two 
sources of uncertainties studied, which was found to be maximally of the order of 10%. Now that, 
thanks to the recent developments in event-by-event hydrodynamics and high-precision data, more 
precise comparisons of flow coefficients between the data and the theory are becoming possible, 
it is important to quantify all the relevant uncertainties to this precision, so that the QCD matter 
properties could eventually be determined from the measured particle spectra and their azimuthal 
asymmetries. 
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